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I. INTRODUCTION 

Brownian motors have been studied intensively since 
the early 1990s [Tl H El H]. This interest coincided 
with developments in bioengineering and nanotechno- 
logy, where understanding and designing a motor in the 
shape of a small biological or artificial device is an im- 
portant issue. Most of the motors investigated in this 
context are powered by chemical energy. Brownian mo- 
tors driven by a temperature gradient [H [H [3 13 [H [10] 
have a fundamental appeal, since their operation is di- 
rectly related to basic questions such as Carnot efficiency, 
Maxwell demons and the foundations of statistical me- 
chanics and thermodynamics [TTl [T2j [131 [HI HI]- The 
additional significance of the thermal Brownian motor 
comes from the recent observation that it can operate as 
a refrigerator [T6l H^, see also [T8^. In fact, this prop- 
erty is, at least in the regime of linear response, a direct 
consequence of Onsager symmetry [19] : if a temperature 
gradient generates motion, an applied force will generate 
a heat flux. This principle is well known in its appli- 
cation to electro-thermal devices, displaying the Peltier, 
Seebeck and Thompson effects [20] . At variance however 
with these macroscopic devices, rectification of nonequi- 
librium thermal fluctuations provide the driving mecha- 
nism for Brownian refrigeration. The latter become more 
prominent, and so do the resulting motor and cooling 
functions, as the apparatus becomes smaller. 

Since the properties of the Brownian heat pump follow 
by Onsager symmetry from those of the Brownian motor, 
we ffist focus on the latter. 



II. BROWNIAN MOTORS 

In earlier suggestions, Brownian motors move linearly, 
which obviously poses difficulties when comparing with 
real systems, or suggesting a technological implementa- 
tion of a Brownian motor. In this paper we introduce a 
Brownian motor, driven by thermal fluctuations, that is 
free to rotate around a fixed axis. Rotational motion typ- 
ically encounters less friction than purely translational 
movement and a rotating force is easier to apply than a 
linear force. We propose a device that exploits the ran- 
dom nature of the perturbations from its environment 
maximally to produce a net directed motion. 

Molecular motors operating within biological cells, al- 
though chemically driven, are also subject to random 
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FIG. 1: Fi-ATPase is a naturally occurring rotating molecular 
motor. It works as a pump for ions through a membrane. 
Its driving force is chemical, through the hydrolysis of ATP. 
Our interest in this molecular motor lies in its environment 
(in the cytoplasm, subject to fluctuations), its construction 
(rotational, and through a biological membrane), its physical 
characteristics (size of the order of 10 nm) and dynamical 
properties (rotational frequency of the order of 100 Hz). 

motion. Fi-ATPase is a well-known rotating motor (see 
Fig. [l] for a diagram). A direct observation of its rota- 
tion, driven by the hydrolysis of adenosine triphosphate 
(ATP), was ffist reported in [2T1[22. Later experiments 
[23] revealed the direction in which the ATP motor spins. 
It is about 10 nm in size and typically rotates with a fre- 
quency of 100 Hz. The observed rotary torque reaches 
more than 40 pN nm. The relation between the geometry 
of the rotating object, specifically its chirality, and its ki- 
netic properties, such as the average motion and friction 
can be of interest to microbiology. One might also imag- 
ine artificial devices inspired by the existing biological 
examples. Proteins could be used as the building blocks 
of mechanical devices and artificial biological membranes 
as means to separate reservoirs and keep them at different 
temperature. Small moving parts in the area of micro- 
electronics are also subject to random fluctuations. 

As a side note we remark that the rotational three- 
dimensional model we will present, can be used to de- 
scribe the essence of the device R. Feynman presented in 
his Lectures on Physics [24 [for a sketch of the ratchet 
and pawl mechanism, see Fig. [2|a)] to illustrate the im- 




FIG. 2: (a) The ratchet and pawl mechanism used by R. Feynman to illustrate the impossibility to extract work from a system 
in equilibrium. An axle with vanes in it (in the right reservoir ) is bombarded by gas molecules at temperature Ti. The pawl 
in the left reservoir (surrounded by a gas at temperature T2) seems to allow only one rotation sense of the ratchet, that is 
connected with the axle. The fluctuations of the paddle in the right reservoir would then be rectified. Via a comparison of the 
probability to move forward and backward, Feynman showed in his Lectures on Physics [24] that at temperature equilibrium, 
Ti = T2, no average motion occurs and the device cannot be used to do work, such as to lift a weight. For Ti > T2 however, 
average motion does take place and the ratchet works as an engine. Feynman noted that for Ti < T2, the ratchet goes backward. 
Note that the rectification manifested by the device originates from the asymmetry of the ratchet and pawl mechanism. Our 
model remains close to Feynman's system, as it can be applied to rotating three-dimensional objects of any shape, while it 
simplifies the asymmetry requirements to the geometrical properties of the device, (b) A rotating Brownian motor where two 
parts in isolated thermal reservoirs (temperatures Ti and T2, particle densities pi and P2) are connected through the axis of 
rotation. Collisions with particles in the reservoirs will cause fiuctuating rotational movement of the heavier motor, which 
under appropriate conditions, will propel the motor with a nonzero average angular velocity. 



possibility of a Maxwell Demon, that would be able to 
extract work from a system in equilibrium. It was also 
shown that the same device can lift a weight (do work) 
at temperature disequilibrium. 

These reasons motivate us to study a chiral Brownian 
motor in detail. The constituting parts of the motor are 
in different thermal reservoirs. In an analytical analysis 
we will derive expressions for the kinetic properties of the 
motor as a function of the external parameters of the sys- 
tem. It will be made clear that temperature equilibrium 
between the different reservoirs results in zero average 
motion and hence prohibits the creation of a Maxwell 
Demon. On the other hand we will demonstrate the im- 
portance of the configuration of the building blocks of the 
motor and their actual shape, as well as their position 
relative to the rotation axis. Some emphasis is put on 
finding optimum operation, yielding maximum average 
angular velocity. The chiral Brownian motor presented 
in this paper is a precursor for the chiral Brownian re- 
frigerator presented in the next. The exact relationship 
derived here between the angular velocity and the tem- 
perature gradient will be a crucial step in the investiga- 
tion of the cooling potential. 

The concrete model we propose consists of at least two 
parts, each residing in a thermal reservoir i = 1,2,..., 
that are rigidly connected with each other through a ro- 
tation axis. Fig.[2|b) shows the construction of two parts 
in reservoirs of temperature Ti and T2, and particle den- 



sities pi and p2- We expect that the fluctuations from 
collisions with particles in the thermal reservoirs will un- 
der certain conditions be rectifled, resulting is an average 
rotational motion, clockwise or counterclockwise. As we 
will show, these conditions are (1) thermal disequilib- 
rium, Ti 7^ T2, and (2) asymmetry (or chirality) through 
the geometrical shape of the motor parts. 



III. FROM FLUCTUATIONS TO THE 
ANGULAR VELOCITY 

Our analysis is based on an exact calculation of the 
probability for the motor to change its rotating speed by 
a certain amount when subject to thermal fluctuations. 
We will show that an exact solution can be reached when 
the fluctuations are in the form of collisions of particles 
of an ideal gas at temperature equilibrium with the sur- 
face of the motor. A master equation for the probability 
density to observe an angular velocity P{u;^ t) at a cer- 
tain time t can be proposed if the particles are presumed 
to collide not more than once, and only with the motor. 
This condition implies that the gas is in the high Knud- 
sen number regime and that the shape of the motor is 
such that recollisions are impossible. We therefore limit 
the parts of the motor to convex and closed shapes. 

We are interested in the case where the motor, with to- 
tal mass M, has no translational degree of freedom and 



3 



a single rotational degree of freedom. Parts of the motor 
reside in different thermal reservoirs but are considered 
rigidly linked. Choosing the z-axis as the axis of rota- 
tion, we can write for the angular velocity = (0,0,a;). 
The inertial moment of the motor with respect to the 
rotation axis is simply denoted as /. 

Under these conditions the probability density P{uo^t) 
obeys a master equation, 

dtP{uj,t) = J du' [W{uj\uj')P{u',t) - W{uj'\u)P{u,t)] , 

(1) 

where W{uo\uo') is the transition probability per unit time 
for the motor to change its angular velocity from uo' to cj. 
The solution is based on the van Kampen l/l^-method 
[25] . A Taylor expansion of the first term of the integrand 
in the angular velocity change, v = uo — uo' ^ leads to 

d,P{u;,t) = Y,^-^i^) {an{u>)P{u,,t)}. (2) 

n=l * ^ ^ 

In this expression the so-called jump moments appear, 
given by 

an{uo) = J v''W{uo;v)dv. (3) 

A notation W{uo';v) = W{uo\uo') is used. With the time 
evolution of the probability density known (Eq.|2|, it is 
possible to derive a coupled set of equations for the mo- 
ments of the angular velocity {uo'^): 

a,K) = f^Q(.;"-'=afcH), (4) 

with (2) the binomial coefficients. Our strategy is now 
clear: first find an expression for the transition proba- 
bility W{uo\uo')^ then calculate the jump moments an{uo)^ 
and finally the moments of the angular velocity {uo'^). 

However, the coupled set of equations Eq. [4] cannot be 
solved unless we expand each equation into powers of 
a small variable, and ignore terms after a certain order. 
For the expansion variable we will use e = tq ^Jrnfl^ with 
ro = \/TJM the radius of gyration of the motor. We also 
introduce an effective temperature Teff, so that to first 
significant order, in the regime of stationary motion, the 
average kinetic energy of the motor is given by 

= ifcsTeff. (5) 

In the calculation it is convenient to do a transforma- 
tion to dimensionless variables, by scaling the angular 
velocity uj and the jump moments as follows: 

An(0 = (6) 

Our selfconsistent definition of the effective temperature 
Teff then leads to (^^) = 1 for the stationary state to first 




FIG. 3: In each reservoir, the motor part can be described by 
the boundary f(x,y), where the rotation axis is at the origin 
of the reference frame. It is convenient to also know explicitly 
the orientation of the boundary at any point. This is given 
by the polar angle cp of the normal outward vector e± on the 
surface at this point. The unit vector is tangential to the 
boundary. 



order in e. The set of coupled equations for the moments 
^^n^ = / ^^P(^,t)d^ remains 

9t{n=p^l^^{C-'A,{0). (7) 



IV. TWO-DIMENSIONAL MODEL OF THE 
MOTOR 

The motor consists of parts with hard surfaces of arbi- 
trary (but convex) shape, each described by their bound- 
ary and inner mass distribution. Many of the important 
features already appear in a simpler two-dimensional sys- 
tem, which we present first. Here the motor consists of 
two-dimensional shapes, each in two-dimensional reser- 
voirs. We choose a cartesian coordinate system as fol- 
lows: the z-axis coincides with the rotation axis, while 
the x?/-plane is parallel to the reservoirs. In each reservoir 
i, the shape of the motor (part) is defined by its bound- 
ary fi{x^ y), given as a vector with the rotation axis as its 
origin (see Fig.[3|. The perimeter of the boundary is de- 
noted Li. Henceforth we will just write r(x, y) for r^(x, y) 
as no confusion can arise in subsequent expressions. It 
is convenient in the further derivation to explicitly know 
the inclination of the boundary at f{x,y), for which we 
use the orientation of the normal outward unit vector 
e_\_ = (cos (p, sin (p) on the boundary, determined by the 
polar coordinate cp from the x-axis. 



A. Conservation rules 

Fluctuations of the angular velocity uo of motor arise 
from collisions on the surface with gas particles of mass 
m. Such a collision - presumed instantaneous and per- 
fectly elastic - changes the velocity of the gas particle 
= {v'^^Vy) into V = (vx^Vy) after the collision, while 
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the motor changes angular speed from co' to co. The iner- 
tial moment of the motor about the rotation axis (z-axis) 
is denoted /, while its mass is M. Conservation of the 
total energy requires that 

^luj''^ ^^mv^'^ ^^mvy'^ = ^lu'^^^mvl^^mvl, (8a) 

while conservation of the angular momentum in the z- 
direction yields 

m{xVy — yv'^) + lix)' = mixvy — yvx) + loo. (8b) 

Also we suppose the interaction force is short-ranged and 
central, implying that the tangent component of the mo- 
mentum of the gas particle on the boundary is conserved. 
Choosing the tangent unit vector e^, = (— ^hup^ cos (/?), so 
that (e^, e^i) forms a positive orthonormal base, we write 

v' ■ e^^ = V ■ e^|. (8c) 

The conservation laws (Eqj8| produce a solution for the 
postcollisional angular velocity cj, 

, 2{uj'y + O cos (/; - 2{uj'x - v' ) sin 
X sm (f — y cos + ^ (x sm (p — y cos (p) 

Introducing 

r,! = r- e^i = —X sin + cos (10) 

and the precollisional speed of the boundary at position 
f, 

V' = u' X r = {-uj'y.uj'x), (11) 
so that we can write 

V'^ = V' ■ = —uo'y cos (p + uo'x sin (/:?, (12) 

and 

/Sy[_ = {V — v')-e_\_ = — {uj'y-\-v'^) cos (p-\- {uo'x — Vy) sin (p^ 

(13) 

the transition in co can also be written as 



(14) 



B. Transition probability 

Next, we set out to find the crucial transition probabil- 
ity W{u;\u;') for the motor to change its angular velocity 
from co' to a; in a unit of time. Every reservoir i contains 
a gas with particle density pi and velocity distribution 
The contribution dWi to the total transition probability 
W{u;\u;') from all possible collisions of particles in gas i in 
a time interval dt on a boundary section of length dli , at 



position r(x, y) and with orientation can be expressed 



as 



xS 



-OO J — oo 

X H[iV' - V') ■ e^]\{V' - V') ■ e^lpiMv'., v'y) 
, 2{uj' y -\- v'^) COS (p — 2{uj' X — Vy) sin (p 



xsinif — y cos (f -\- — {x sin (f — y cos (f) ^ 



(15) 



with H the Heaviside step function and S Dirac's dis- 
tribution. We multiplied the particle density pi with the 
volume of the gas that is passed by the boundary element 
dli in a time unit, considering only those gas particles 
that comply with the collision rules. This can be written 
in short form as 



dv'^ / dv'yH[/\Vl]\/W^ 

-oo J — oo 



X pi(j)i{v'^,v'y)5 



uo — uo 



■ (16) 



The total transition probability is then found by inte- 
grating over all boundary elements dli and summing over 
all reservoirs: 

W{u\uj') = ^(f dWi{u;\u;'). (17) 

^ ^boundary 

Henceforth we will simply write f when we imply the line 
integral over all boundary elements. 

For a Maxwellian velocity distribution at temperature 
Ti. 



(t>i{Vx,Vy) = 



2iTkBTi 



exp 



-m{vl + vl 
2kBTi 



(18) 



the integrals over the speed of the colliding particles can 
be performed explicitly, resulting in 



2TTkBTi 



X exp 



X (r„ + fHUuj - uj')rJ\uj - uj'l 

m f, I .oj' — u! I 



• (19) 



C. Moments of the angular velocity 

Now that we have obtained an exact expression for the 
transition probability W{uo\uo')^ we turn our attention to 
the jump moments. 



W(cj; v)dv, 



(20) 
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and then the moments of the angular velocity. Careful 
consideration of the sign of r,, + I /mr^^ in 



m 



X H 



[v] I dk^H[-v] [ dk] \v 



X exp 



, (21) 



where v = a; — cc;' is the change in angular veloc- 
ity, leads to an exact expression for the jump mo- 
ments. In terms of parabolic cylinder functions, D^(z) = 
(exp[-z^/4]/r[— n]) exp[-2:x - x'^ /2]x~'^~'^dx (for 
n < 0) the results are 



an{uj) 



X (b dliPi r,| + 



m 
I 



X exp 



^2, ,2 



D-n-2 



(22) 



Rescaling the jump moments using dimensionless vari- 
ables ^ = Lo^jl /kBTes and e = tq ^Jm/I, where Tq = 
I/M, leads to 



2^ 
X (h dl 



X exp 



m vTeff 



i/2 



2^2 





"n + 2" 




2 



2 y.-(n,/ro)^e 
n + 2 1 e^T;^ 



(r„/ro)^C^ 



V Ti To 



' 2 ' 2 

"n + 3 



X $ 



n + 3 3^ £2 Teff 

2 ' 2' 



(r„/ro)'e' 



(23) 



Here ^ represents Kummer's function [26 . 

We can express both the exponential function and 
Kummer's function in a power series, 

v2 ^3 



exp[z] = 1 ■ 



z z 
1! ^ 2! 



z 
3! 



(24) 



^a;j;z] = iFi[a;j;z] = 1 



a z a{a + 1) 



7I! 7(7 + 1)2! 
a{a -h l){a + 2) z^ 



7(7 + l)(7 + 2) 3! 



, (25) 



Considering that the parameter e = ^Jm|M is small for 
gas particles much lighter than the motor, we arrive at 
a series expansion for the jump moments in e. We sub- 
stitute this expansion in the set of equations (Eq. 23 ) 



coupling the jump moments A^(0 with the moments of 
the angular velocity (^^). For n = 1, with r = e^t, this 
results in 




dk 



6(0 



T 

-J- 1 

^eff 



7, 



ksTj 
m 



V TT 



3 



dk 



dli 



The term in £ ^ disappears because 

dliV^^ = ^ dli ■ f = J (V X r) • e^dAi = 0. 
Similarly for n = 2, 



(26) 
(27) 



12 



m 



dli 



/^eff . 



-teff V J-i 



dk I - 



T 

TT V Teff 

Ti~ 



He) 



3 Ti^^ ^ 



dli 



dl 



(28) 



To lowest order in s we can extract from Eq. (26) a 
linear relaxation law for rotational movement, IdtJuj) = 
r/, describing a net frictional torque r/ exerted on the 
motor as a result of all collisions .With r/ = —j{uj) and 
7 = li we derive a microscopic expression for the 
friction coefficient 7^ of each part of the object: 



7i 



kBTifn 
27r 



(29) 



To order the average angular velocity in a stationary 
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FIG. 4: In each reservoir, the motor part can be described by 
the boundary f{x,y), where the rotation axis is at the origin 
of the reference frame. It is convenient to also know exphcitly 
the orientation of the boundary at any point. This is given 
by the polar angle cp of the normal outward vector e± on the 
surface at this point. The unit vector el, is tangential to the 
boundary. 



e± = (sin6>cos(/?, sin6>sin(/?, cos6>), described by two an- 
gles and polar and azimuthal angles in a spherical 
coordinate system. 

The assumption that there is only a central force dur- 
ing an interaction with a gas particle, entails that there 
is conservation of momentum of a gas particle along any 
tangential direction. So for the speed of the gas particle 
before ( v' = {v^^ , Vy ^ v'^)) and after {v = (vx^Vy^Vz)) the 
collision, we can write this condition formally as 



V • e/c,, 



V • e/c,, 



1,2, 



(33) 



where ei^,, and 62,11 are two distinct unit vec- 
tors perpendicular to e±. It is convenient 
to use ei^ii = (— sin(/:), cos(/:?, 0) and 62,11 = 
(cos 6 cos (p^ cos 6 sin — sin 6) . Together with con- 
servation of total energy and angular momentum in 
the z-direction (the expressions are the same as in the 
two-dimensional analysis) we find a relation for the 
change of angular velocity induced by one collision: 



state appears from Eq. (26) as 



4/ EiPiVTi§dk 



kBTn 



dhrt 



with second moment, 



B-Le« 



(30) 



(31) 



while the effective temperature is found according to its 
definition, 



Teff=(5^7.T,)/(5^7.)- 



(32) 



Using higher order terms in the expansions for (^^) re- 
sults in correction terms to the expressions for (uj) and 
(cj^). The second terms are in both a factor m/M smaller 
than the first terms. 



V. THREE-DIMENSIONAL MODEL OF THE 
MOTOR 

The results of a fully three-dimensional analysis are 
very similar to those derived in the previous two- 
dimensional case. We will clarify the key differences here. 

The motor parts are now determined by their surface 
Si {i is the reservoir in which the part resides), to be de- 
scribed by a vector r(x, z) from the origin. Again we 
choose the rotation axis to be along the the z-axis. The 
orientation of the surface at a certain location (x, z) is 
uniquely determined by the normal outward unit vector. 



I 



(34) 



{V^ — v^) ■ e± is the component of the velocity difference 
between motor and gas particle perpendicular to the sur- 
face at the place of impact, r,, is now defined as 



-X sm 6/ sm (f -\- y sm tl cos (p = sm tl r • ei, 



(35) 



r,! is zero in locations where the surface is perpendicular 
to r, these coincide with zero momentum transfer. Max- 
imal momentum transfer and r,, occurs when the tan- 
gential plane to the surface at this location crosses the 
rotation axis. 

The transition probability dW{u;\u;') caused by all pos- 
sible collisions with a surface element dSi of the motor, 
is then found by integrating over all velocities that obey 
the collision rules. 



dvl 



/-j-CO p-rOO n- 

dv'^ / dv'y / 
-00 J —00 J — 

X H[{V' - v') ■ e±]\{V' - v') ■ e±\ 
Ai ' ' ' o(^'-«'*')-e± 



r„ + 



(36) 



Adding the contributions of all surface elements dSi in 
all reservoirs i, gives us the total transition probability, 

W{uj\uj') = ^ [ dWi{u;\u;'). (37) 

^ J surface 

Again, for a Maxwellian velocity distribution. 



m 



27rkBTi 



3/2 



exp 



2kBTi 



(38) 
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we can do the integration over v^jVy, analytically and 
find 



B. Chirality 



W{io\u;') = ^ E / "^^iPi 



27TkBTi 



X (r„ + fH[{co - U!')r„]\u! - co'\ 

mr,. 



X exp 



m 



2kBTi \^''"^ mr?' 



(39) 



This expression is identical to its two-dimensional equiv- 
alent (Eq.[l9|, apart from the different definition of r,,, 
and obviously an integration over the surface instead of 
the boundary. The previous algebraic technique can then 
be applied to derive results for a general shape of the mo- 
tor, such as for the average angular velocity in a steady 
state, 



^/2'!rkBm pi{TeS - Ti) J dSi\ 



and the friction coefficient, 



27r 



(40) 



(41) 



(42) 



where Teff is still defined as 

i i 

VI. ANALYSIS AND DISCUSSION 



Now that we derived analytical results for any shape 
and any number of reservoirs, we are ready to analyze 
concrete systems. We are interested in the role of ex- 
ternal parameters, such as the temperature and the den- 
sity of the gas, and in the construction and shape of the 
motor itself. Much of the analysis can be applied to the 
simpler two-dimensional case, but references to the three- 
dimensional case are made where they are appropriate. 



The next element 
fdkr^^ in Eq. [43 
reservoir i that is 



we want to discuss is the factor 
Consider a motor shape in one 
symmetrical with respect to a plane 
through the rotation axis. A simple argument reveals 
that § dlivf^ = 0: for every point (x, y) on the boundary 
of the shape with value r,, there can be found a point 
{x'^y') for which r,, = — r,,. The contour integral of r,, 
is therefore zero, considering that the line element dli is 
positive. 

A construction that consists entirely of symmetric 
shapes will yield zero average rotation. Such a construc- 
tion in its most simple form could consist of flat blades 
through the rotation axis in every reservoir. To find a 
net angular velocity, the motor must have at least one 
chiral part. The factor f dlirf^ will be analyzed in more 
detail in a later section, and we will show that under cer- 
tain conditions it can also become zero even for a chiral 
configuration. 



C. Friction and propulsion 



In the full expression 



(44) 



the factor § dlirf^ in the denominator stems from the fric- 
tion each motor part encounters while rotating in the gas. 
If we look at optimizing the motor, the first idea would 
be to minimize this factor. A surface where r,, is zero at 
every point corresponds to a sphere, but ^ dkrf^ will be 
zero as well, resulting in zero net motion. Large average 
angular velocities will be obtained then by a compromise 
between a small ^ dlivf^^ and a large § dlirf^. The propul- 
sion of the motor originates in the factor f dlir^. The 
largest friction will be experienced by shapes where r,, is 
maximal. This corresponds to a (flat) surface, or blade, 
through the rotation axis. 



A. Temperature gradient 

When the thermal reservoirs are at equilibrium with 
each other, we immediately see from Eq. 32 that Ti = 
T2 = • • • = Ti = Teff, independent of the construction we 
propose. The average angular velocity 

kBTTl 



^Pi{T,f,-T,) ^dhrl (43) 



becomes zero. It is impossible to extract net motion from 
a system in equilibrium. 



D. Motor configurations 

We turn to the question of how to configure the motor. 
Leaving the exact choice of the shape for later, we tackle 
the following the question: if we have a certain part of 
the motor in one reservoir, how will the placement of 
the part in the other reservoir effect the motion of the 
motor We start by proposing three simple constructions 
(see Fig.[5|: 

1. The shapes are identical in both reservoirs 
(Fig.|5^). This includes the location of the rotation 
axis with respect to the shape. The exact shape 
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FIG. 5: Three possible configurations of the motor parts in 
a two-reservoir system. In (a) the part in the first reservoir 
is copied exactly to the other reservoir. Not only the shapes 
are identical, also their position relative to the rotation axis. 
In (b) the motor part in one reservoir is reflected in the other 
reservoir. The reflecting plane passes through the rotation 
axis. In (c) a general shape is combined with a blade: a plane 
of length L (and height H in three dimensions). Note that the 
system is rotationally invariant in each reservoir separately. 



can be determined afterwards. With 

dhr^^ = ^dl2r^= jdlrl, (45) 
dhrl = (f dl2r^,= l dlr^, (46) 



Eq. 44 simplifies to 



PiP2(T2 - tI^\T2 - Ti) § dl 



4/ 



§dlrY 
(47) 



2. The shapes are still general but they are exact mir- 
ror images of each other in both reservoirs (Fig.jsjD). 
Also the location of the rotation axis with respect 
to the shape is mirrored. The mirror axis (plane) 
is through the rotation axis, but its orientation is 
of no importance, as our system is rotationally in- 
variant. Writing 



dlir'f^ 
dhrl 



dhr'^ 



dUr^r- 



dlr'^^^ 
dlrl, 



(48) 
(49) 



Eq. 44 now becomes 



4/ {piTl'^ + P2TI'^Y §dlrl' 

(50) 



3. We use a general (yet unknown) shape in the first 
reservoir, while in the second reservoir we put a 



blade of length L, rotating about one end (Fig. [5]: 
Omitting the index i = 1, and identifying 



j dl2rl = 2LV3, 
dhr^, = 0, 



(51) 
(52) 




FIG. 6: A motor with multiple identical shapes in each reser- 
voir. For the depicted configuration with equal multiples, 
n = m = 4, we show that the resulting average angular veloc- 
ity is the same as that for only one shape in each reservoir. 



we obtain 



4/ (piT//' / dl rl + p2T2'/'(2L3/3))2 ' 



(53) 



Comparing the three suggested configurations, we see 
configuration (1) is even when the temperature differ- 
ence AT — T\—T2 is inverted, while (2) and (3) are odd. 
For small temperature differences, {uj) is approximately 
parabolic in AT, while (2) and (3) are linear. For small 
AT therefore (1) yields much lower angular speeds than 
(2) and (3). 

For a small temperature difference a rather technical 
calculation shows that the average angular velocity for 
configuration (2) is at least twice that of construction (3) 
for the same general shape with similar linear dimensions 
as the blade. 

In a numerical procedure described later, we discovered 
that the configuration of two mirror shapes in the two 
reservoirs produces the maximal average angular velocity. 
In other words, given a certain part in one reservoir, the 
highest average angular velocity is obtained by using the 
reflected shape in the other reservoir. 

One could think of multiple but identical structures 
(blades) in each reservoir, as illustrated in Fig.|6] Ignor- 
ing the increased probability of multiple collisions of gas 
particles with the motor, our theory leads to an average 
angular velocity 
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^/2^^ P,P2{T2 - Ti) (t^/^ § dh rl § dh rl - tI''' § dh § dh r\ 
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(54) 



for a system with n identical blades in reservoir 1 and 
m identical blades in reservoir 2. The appearing contour 
integrals are over one shape of the set of identical shapes. 
For the same number of shapes in both reservoirs, n = m, 
the average angular velocity is the same as with only one 
blade in each reservoir. The result also shows that it is 
beneficial to have the highest number of blades in the 
reservoir with the highest pT^/^ § dl rf^ factor, or simply 
the highest pT^/^ factor if the blades have the same shape 
in both reservoirs. 



E. Globular proteins 

Looking for real- world candidates to fill the role of our 
Brownian motor, we turn our attention to biological sys- 
tems. In the further analysis we want to use physical 
values for the dimensions, masses and so on. A possi- 
bility is to apply our model to globular proteins, which 
could give the shape of the motor parts. The two parts 
would reside in a water environment, separated by a lipid 
membrane. 

To obtain orders of magnitude for our results we will 
refer often to the values in TableU 



Mass of one part M/2 1.66 X 10~^^kg 

Density of the motor pm 1380 kg m"*^ 

Volume of one part V 120 nm^ 

Radius of one part (if assumed spherical) R 3 nm 

Particle mass (H2O) m 2.992 x 10~^^kg 

Reservoir temperature Ti , T2 ±300 K 

Reservoir particle density Pi?P2 ±3.3 X 10^ 




FIG. 7: The temperature dependency of ooq for a configuration 
of (a) identical shapes and of (b) mirror shapes in the two 
reservoirs. One reservoir is kept on a fixed temperature Ti = 
300K, while the temperature T2 in the other reservoir changes 
in a range of 300 =b IK. For (a) we see a nearly parabolic 
dependency on T2 — Ti. The sense of rotation remains the 
same if we switch part 1 and part 2 of the motor, and the 
angular velocity is rather small. For (b) we see a much larger 
effect, a nearly linear dependency on T2 — Ti, and the sense 
of rotation is inverted by switching part 1 and part 2 of the 
motor. The densities are taken the same in both reservoirs. 



motor part, is multiplied to make S scale-invariant. S 
will be discussed in detail in the next section. Eqs. 47 



50 can then be written as 



(56) 



n28^ 



What remains is a factor loq that depends on the specific 
configuration, the reservoirs temperatures and densities, 
and the masses of the motor and particles: 



TABLE I: Typical parameters used for the motor and en- 
vironment. The values of the individual parts correspond to 
those of globular proteins. 



V2^k^ PiP2{Tl^^ ^tI^\T2 



4MVA 



1/2 



(57) 



F. External parameters 

In the cases of identical shapes and mirror shapes 
(Eqs. 47 , 50 ) we can separate from the expressions for 
{co) a shape-dependent factor. 



TJdfi 



.2 • 



(55) 



A factor 
mensions) 



A (or V^/^ if we prefer to work in three di- 
which stands for the typical dimensions of a 



We have used the notation uJq for the configuration with 
two identical shapes in the two reservoirs and Uq for the 
configuration where the shapes are mirror images. 

coq is also dependent of the size of the motor. Because 
the mass M of the motor is also size-dependent, the full 
dependency could be written as M\fA = 2p^A^/^ in two 
dimensions, or MV'^^^ = 2pmV^^^ in three dimensions, 
if the density pm of the motor interior is considered con- 
stant. Therefore ujq oc M~^/^ in two dimensions and 
coq (X M""^/^ in three dimensions. 

Fig. [T] shows the temperature dependency of cjq, cal- 
culated with the physical values of Table [l| and equal 



reservoir densities. As mentioned before cOq 
T2 — Ti, while uJq is quadratic. 



is linear in 
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FIG. 8: Three simple shapes for the motor parts. All all 
determined by a minimal number of parameters. The shape 
is determined by the angle a, the size by one side R, and 
the position with respect to the rotation axis (at the origin) 
by one coordinate {x,y). For the triangles (Motor 1: right 
triangle; Motor 2: isosceles triangle) (x^y) is the location of 
the center of mass, while for the disk sector (Motor 3) (x, y) 
corresponds to the center of the disk. 



G. Shape factor 

Next, we consider the size-independent geometrical 
factor, 



(58) 



which is comprised of an interior factor M//, and an 
exterior (boundary) factor § dl rf^ / § dl r,^ 

The factor M/I is actually independent of the mass of 
the motor because the inertial moment / is proportional 
to the mass M. It only describes the spacial distribution 
of mass. For a homogeneous motor interior it is given by 



M _ 



(59) 



The integral is over the entire interior of the motor, and 
r is the distance of an interior point to the rotation axis. 

Finally the factor § dlrfj § dl rf^ depends on the exact 
form of the boundary of the motor parts, where r,, is to 
be measured from the location of the rotation axis. The 
integrals are over the entire boundary. 

To enable us to get an understanding of the geometri- 
cal factor, we introduce three simple realizations (Fig.[8| 
(in two dimensions and with homogeneous mass distri- 
butions), of which the boundary can easily be described 
analytically. For these prototype shapes all factors can be 
expressed in closed form. The shape of the motors parts 
are respectively a right triangle (Motor 1), an isosceles 
triangle (Motor 2) and a disk sector (Motor 3). Both 
the dimensions of the motor parts {R) and the shape 
(a) are fixed with one parameter, making a comparative 
study easier. We would also like to specify the location 
of the motor part relative to the rotation axis with one 
representative point with coordinate {x,y). For motors 1 
and 2 we choose the center of mass and for motor 3 the 



Motor 1 
Motor 2 
Motor 3 



I/M = R^/18 + x^ + y^ 

I/M = R^{2- cos 2a) /18 + x^ 



I/M : 



2 



4:Ry sin a 
3^ 



TABLE II: The ratio of the inertial moment / over the mass 
M of the motors in Fig. [S] given as a function of the location 
{x,y) (of the center of mass for Motors 1 and 2, and of the 
center of the disk for Motor 3) with respect to the rotation 
axis and the shape parameters, angle a and size R. The dis- 
tribution of mass within the motors is assumed homogeneous. 



center of the disk sector (see (Fig.[8|) as this represen- 
tative point . For these simple motors we can calculate 
analytical expressions for §dlrf^ and §dlrf^. As these 
expressions are rather long, and the details are not of 
immediate relevance, we have put them in the appendix. 
As an example consider the Motor 2 case: 



2R f 

dl r,^ = (R^^ 3Ry cos 3a + 9y^ + 9x^ sin a 
sin^ a [R^ (3 sin a - 2 cos 2a) + 9x^-9^^) ), (60) 



and 



R 

dlr^^ = —X sin 2a ( R{3y + R cos a) (1 — 2 cos 2a) 



+ 3 (x^ - 3^^) cosaj. (61) 

The intertial moments / for the three motors are given 
in Table [m Some physical properties are immediately ap- 
parent from these expressions. For example § dlrf^^ and 
hence the angular velocity (cj), is zero when 

• X = 0: this is when the rotation axis is on the 
symmetry axis of the motor; there is no preferred 
sense of rotation, 

• sin 2a = 0ora = 0ora = 7r/2: the motor is bar 
shaped, and loses its asymmetry (or chirality). 

Note that f dlrf^^ which also appears as a factor in the 
expression for the friction coefficient, is not zero if the 
shape is bar shaped (or symmtrical in general). 

In general the §dlrf^ factor (Eq.pl] for Motor 2) de- 



scribes the asymmetry of the motor^Lt also determines 
the sense of rotation. For example for Motor 2, the rota- 
tion sense is inverted when the rotation axis is placed on 
opposite sides of the symmetry axis, x = xq and x = —xq. 

More features can be seen from plots of the angular 
velocity of the motors as a function of their shape and 
configuration, see Figs.|9] (Motor 1), [To] (Motor 2) and 
11 (Motor 3). Note that a coordinate change was made, 
X —X, y —y. This means the coordinate (x, y) in the 
plot corresponds to the location of the rotation axis with 




FIG. 9: Contour plots of the average angular velocity (oo) (x, y) 
of Motor 1 as a function of the location of the rotation axis 
for several values of the angle a. The center of mass 
of the motor is in the origin (0,0). Regions in black (and 
white) correspond to locations for the rotation axis that yield 
the highest (co) (but in opposite sense). If the rotation axis 
is put on a red curve there is zero average rotation. Maxima 
in {(jj){x,y) are marked by dots, the purple dot reveals the 
optimal place for the rotation axis (the one that gives the 
highest angular velocity). Note that one red curve for the 
a — 45° realisation corresponds with a symmetry axis of the 
shape. The motor will not show directed motion if the axis 
is place there. Note also that mirror shapes (such as a = 20° 
and a = 70°) show opposite rotation sense, for all locations 



respect to the representative point of the motor (center 
of mass or center of the disk) , which is put in the origin of 
the plots. Figs.|9j [Toj [lT]of {uj){x^y) show lines of equal 
angular velocity (in black) and lines of zero average angu- 
lar velocity (in red). Highest angular velocities are found 
in the black and white regions (but with opposite rota- 
tion sense). Local extrema of {uo){x^y) are represented 
by a green dot while a purple dot is the optimal location 
of the rotation axis. The shape of the motor is drawn in 
yellow. 

We see that the red curves that signify zero average 
rotation can be straight lines when they correspond to a 
symmetry axis of the shape (the ?/-axis in Figs.jlO, 11 for 
all shapes a, but also in Fig. [9] for a = 45° for example), 



but in general they follow a curved path. The regions 
of opposite rotation sense, separated from each other by 
the red curves, form not so trivial patterns. 

Also interesting to note is that the location of the rota- 
tion axis that yields the highest rotation speed is always 



FIG. 10: Contour plots of the average angular velocity 
{uj){x,y) of Motor 2 as a function of the location of the ro- 
tation axis (x,^), for several values of the angle a. For a 
technical explanation see Fig.|9] The symmetry of the shape 
is reflected in the {oo){x, y) plot, in particular in the locations 
for the rotation axis that correspond to zero average angu- 
lar velocity (red curves): the ^/-axis for all the configurations 
and three symmetry axes for a — 30°. Note that the rotation 
sense is opposite for locations on opposite sides of a symmetry 
axis. 



in the vicinity of the rotating Brownian motor and often 
in its interior. Remember that the geometrical results are 
scale-invariant, and the relative locations of the maxima 
(and zero lines) are independent of the dimensions of the 
motor. 

For certain choices of the shape and especially of the lo- 
cation of the rotation axis the average rotation speed can 
become zero. It is therefore sensible to investigate which 
configurations yield the highest rotation speed. For the 
three simple motor realizations, we determine the loca- 
tion of the rotation axis that yields the highest shape 
factor S for every value of the shape parameter a, as 
shown in Fig. 12 The angle a that results in the highest 
S is listed in Table ITITl for each of the motors. The cor- 

Considering 



responding shapes are depicted in Fig. 13 



the constraints put on the shape. Motor 1 and Motor 2 
adopt very similar configurations, while the best (convex) 
shape for Motor 3 is a semi-disk. 



H. Optimal shape 

The three motor realizations show that the angular 
velocity is sensitive to the precise shape of the motor. 
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FIG. 11: Contour plots of the average angular velocity 
{(jj){x,y) of Motor 3 as a function of the location of the ro- 
tation axis (x,^), for several values of the angle a. For a 
technical explanation see Fig. [9] 
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FIG. 12: For the simple motor types (Fig.[8|, we determined 
the position of the rotation axis that maximizes the average 
rotation speed for every shape, given by angle a. The values 
of a for which the shape factor S reaches a maximum are 
given in Table [llll It becomes clear that the shape is a key 
factor in the operation of the motor. 



We are interested to know what happens if we relax the 
shape constraints while optimizing for maximum rotation 
speed. We solve this problem using a numerical proce- 
dure. 

The boundary of the motor is modeled as piecewise 
linear. It is defined by the location of n vertices. The 




FIG. 13: The choices for the shape and location of the rota- 
tion axis for each of the three motor realizations that result 
in the highest average angular velocity of the motor. The lo- 
cation of the rotation axis is marked by a star. Although the 
initial constraints for Motor 1 and 2 are different (respectively 
the shape of a right and isosceles triangle), in their optimal 
configuration they are very similar. Motor 3 is optimal in the 
shape of a semi-disk, while we excluded the possibility of a 
concave shape. 



Motor 


Angle a 


shape factor S 


Motor 1 


43.2° 


- 0.618 




46.8° 


0.618 


Motor 2 


17.7° 


0.465 




47.7° 


0.627 


Motor 3 


22.4° 


0.638 




90.0° 


1.54 



TABLE III: For each of the three simple motor protypes, we 
shape factor S that corresponds to the optimal settings (shape 
and location of the rotation axis) is listed. The related shapes 
are shown in Fig.|13| 



numerical procedure finds the optimum location of the n 
vertices, yielding maximum angular velocity, under the 
constraints that (1) the mass M remains constant, (2) 
the shape remains convex. The mass constraint for a ho- 
mogeneous mass distribution translates into conservation 
of total area A. For low numbers n = 3,4,5,6 the op- 
timum location of the vertices is shown in Fig. [14] Note 
that the rotation axis is still fixed in the origin (0,0). 

By increasing the number of vertices n, the piecewise 
linear shape approaches the smooth boundary that yields 
the highest angular velocity possible. In Fig.[T5]the shape 
factor S is plot against the number of vertices n. We 
see a convergence for large n. For n = 100 the value 
of S is 2.29. This is a factor 3.65 higher than the best 
value for a triangular shape, n = 3, S = 0.63. The 
corresponding shape (for n = 100) is shown in Fig. 16 
A tentative explanation for the optimum spiral shape is 
that it combines a long curved section with small r,, (and 
hence small friction) with a short section that is linear, 
providing the necessary propulsion. 

The chirality of the shape determines the rotation 
sense. A motor that consists of two identical optimum 



shapes as shown in the large figure of Fig. 16 actually 
has a negative S. This means the motor rotates clock- 
wise (from y-ax.is to x-axis). Its enantiomer (the small 
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FIG. 14: A numerical procedure was used to find optimum 
piecewise linear shapes, by determining the best location of 
each vertex. The results for a small number of vertices, n = 
3, 4, 5, 6 are shown. Note that the results are invariant for 
rotations with respect the rotation axis, which is located in 
the origin (0, 0) of the coordinate system. 
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FIG. 16: The (two-dimensional) shape of one motor part that 
produces the highest average angular velocity, found with a 
numerical procedure. The rotation axis is marked by a star (at 
(0, 0)), while the center of mass is marked by a dot. The shape 
in the larger figure rotates from the ^-axis to the x-axis, while 
its mirror shape (enantiomer) in the smaller figure rotates in 
the opposite sense. The area of the shape is normalized to 1. 
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100 



the same shape factor S (see Eqs. 47 and 50). Then we 
extended the numerical procedure so that the shapes in 
each reservoir could develop independently. For a small 
temperature difference (AT = 1 K) the shapes become 
almost exactly each others mirror image (area difference 
A1/A2 = 1.00045). Even for a large temperature dif- 
ference (AT = 100 K) we find mirror shapes with only 
a small area discrepancy {A1/A2 = 1.017). The larger 
shape is in the reservoir with lower temperature. To 
obtain a positive rotation sense (from x-axis to ?/-axis) 
of the motor for a construction with mirror images, the 
shape with positive S needs to be placed in the cold reser- 
voir , while its mirror shape with negative S (but equal 
in absolute value) should be in the warm reservoir. 



FIG. 15: The shape factor <S is a scale-invariant measure of 
the influence of the precise geometry of the motor. For <S = 0, 
the average angular velocity is zero. A numerical procedure 
was used to obtain the optimum shape (with highest S) for 
a motor shape that is piecewise linear, i. e. consisting of n 
vertices connected by straight lines. For high n we find an 
approximation of the exact optimum shape and a lower limit 
of the maximal S that can be obtained. S converges to a value 
of about 2.29. A negative S is possible, but this corresponds 
to a shape that is the mirror image of the shape with opposite 
S. 



figure) has positive S and rotates counterclockwise (from 
X-axis to ?/-axis). 

We initially applied the numerical procedure to iden- 
tical shapes in the two reservoirs. We knew we would 
simultaneously find the optimum shape for the construc- 
tion with mirrored shapes in both reservoirs as they share 



I. Physical estimates for the Brownian rotor 

We have now collected all the necessary elements to 
estimate the physical properties of our motor. In a real 
world example of course many of the features discussed 
in this paper will only be of qualitative use. 

We learned that is is advantageous to use a configu- 
ration where the shape of the motor in one reservoir is 
the mirror image of the shape in the other reservoir (sec- 
tion jVIDl). We also found the optimal individual shape 



(section |VIH ). In the results we present here we assume 



these optimizations can be approached. 

When we use the physical properties of small protein 
structures in an environment of water molecules (section 
VIE) as the separate building blocks of our motors we 



arrive at an average angular velocity of about 1500 Hz 
when driven by a temperature gradient of 0.1 K (Ti = 
300 K, T2 = 300.1 K). This corresponds to about 230 
rotations per second. 
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VII. CROSS PROCESSES 

We discussed a Brownian motor, and derived a rela- 
tionship between its motion - the average angular ve- 
locity (uj) - and the applied temperature difference AT. 
This relationship is an example of a cross process. The 
normal process that would give rise to a motion {co) orig- 
inates from a mechanical force. In our system with only a 
rotational degree of freedom this mechanical force would 
be in the form of a torque r along the z-axis. 

Cross processes are very common in physics. One well- 
known example is the Seebeck effect [20] , where a temper- 
ature difference over an electric conductor causes an elec- 
tric current. The Seebeck effect has a reverse or mirror 
cross process: the Peltier effect. Here, an applied electric 
current causes a temperature difference. Processes and 
their mirror processes are related through a general prin- 
ciple of stability. In the example of the Seebeck-Peltier 
effects, consider a system that is originally in equilib- 
rium. If it is perturbed by the application of a temper- 
ature gradient, currents will start to flow (the Seebeck 
effect), which in turn will give rise to a counteracting 
temperature difference (the Peltier effect), attempting to 
cancel out the original cause of the disturbance. Another 
example is Lenz' law in electromagnetism. Moving a fer- 
romagnetic core into a coil will induce currents in the 
coil. These currents are such that the resulting magnetic 
field will expel the core, hence counteracting the original 
disturbance. 

Now that we established the existence of a cross pro- 
cess in the Brownian motor system and showed the rela- 
tionship between a cross process and its mirror process, 
the question naturally arises: what is the mirror process 
in our system? It should be one that counteracts the 
original perturbation - a temperature difference between 
the reservoirs. A flow of heat Qi^2 from reservoir 1 to 
reservoir 2 would do exactly that. And we know the nor- 
mal process (that causes rotational motion) is induced 
by a torque r. The mirror cross process then is a heat 
flow caused by a torque. When we perturb our system 
in temperature equilibrium by applying a torque r, the 
motor will of course rotate. How can the system react to 
counterbalance this motion? By creating a temperature 
gradient between the two reservoirs, that according to 
the theory of the Brownian motor will cause rotational 
motion. The sign of the temperature difference is such 
that the resulting rectified Brownian motion opposes the 
motion started by the torque. 




FIG. 17: The Brownian refrigerator is a rotating device con- 
sisting of two connected bodies, each in a temperature reser- 
voir. By applying a torque on the apparatus, a heat flow 
Qi^2 will arise that can cool down reservoir 1 at the expense 
of heating up reservoir 2. It is possible to obtain the condi- 
tions where this cooling power is larger than the dissipative 
heating Qji. 



each reservoir. The parts are rigidly connected via a ro- 
tation axis and are subject to random collisions with par- 
ticles in their reservoirs. These have mass m, considered 
much smaller than the total mass M of the rotating de- 
vice. On applying a torque r along the z-axis the system 
will develop a heat flow (3i^2, cooling reservoir 1 down 
at the expense of heating up reservoir 2. The following 
derivation of this heat flow is done for a construction 
where one part of the device is the reflected copy of the 
other (as in Fig.jsjD) because it shows a linear response 
in small temperature differences Ti — T2 and it produces 
the highest rotating speeds and resulting heat flow. 



IX. LINEAR RESPONSE AND ONSAGER 
SYMMETRY 

Previously we derived a relation between the average 
angular velocity {co) of the rotating brownian motor and 
the temperatures of the two isolated reservoirs, Ti and 
T2, correct to order m/M: 



VIII. MODEL OF THE BROWNIAN 
REFRIGERATOR 

We propose a cooling apparatus or refrigerator based 
on the Brownian motor described before. Two reservoirs 



(see Fig. 17) are thermally isolated from each other and 
initially at temperature equilibrium, Ti = T2 = T. The 
refrigerator consists of a rotating device with a part in 



Tl'\T2-Ti)§dlTl 
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1/2 



(62) 

Here / is the inertial moment of the motor with respect 
to its rotation axis and § dl rf^ / § dl r^ is a geometrical 
factor defined by the shape of the rotating motor parts. 
The geometrical factor is zero in a symmetrical config- 
uration, hence the importance of choosing appropriate 
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asymmetric or chiral elements. For a more detailed dis- 
cussion of the geometrical factor we refer to Section [VI G 



For a small temperature difference AT between the two 
reservoirs, 

Ti = T - AT/2, AT = T2 - Ti < T, (63) 
T2 = T + AT/2, T = (Ti + T2)/2, (64) 

the mechanical response {co) is linear in AT to very good 
approximation: 



V^TTk^rm P1P2 §dlrl AT 
21 (Pi +/02)2 §dlrlT^' 



(65) 



An elegant way to calculate the cooling potential of 
our system is by making use of Onsager's relations [19]. 
We will identify a flow and a force for the two cross pro- 
cesses involved. The alternative is to revisit the analysis 



of Section III , adding an extra torque term to the master 
equation. For the mechanical process we identify a flow 
Ji = {(jj) and a thermodynamic force X2 = AT/T^ in 
the linear relation of Eq.[65] The proportionality con- 
stant T12 of the first Onsager relation. 



is, for our particular system, given by 



^12 



21 



^3/2 . 



(66) 



(67) 



For the second cross process, the cooling process, we 
can identify a heat flow J 2 = Qi^2- The force Xi is given 
by the chemical potential associated with the particle 
flow Ji = (co) of the normal process, and is generated by 
the application of the torque r. More precisely, Xi = 
r/T. Again we expect a linear response 



J2 — T21X1. 



(68) 



Onsager symmetry now tells us that the two propor- 
tionality coefficients of the cross processes are identical, 

T21 =Ti2. (69) 

The heat flow from one reservoir to the other in the linear 
regime now becomes obvious. 



Ql^2 = 



V^TrkBmT pip2 § dl i 
27 {P1+P2Y §dh 



T. 



(70) 



The complete Onsager relations, combining normal 
and cross processes, are given by: 

Ji = LiiXi + T12X2, J2 = L21X1 + L22X2. (71) 

The Onsager coefficients we have not identified yet are: 



Lu = T/7, T22 



7172^^^^ 
7/ ' 



(72) 



while Ji = (cj), J2 = Qi^2, Xi = r/T, and X2 = 
AT/T^ as before. Lu can be associated to the direct 
mechanical response of the motor to the application of a 
torque, while T22 is related to heat conductivity, being 
the coefficient between the heat flow and the temperature 
gradient. 
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FIG. 18: Physical values of the cooling power Qi^2 of the 
Brownian refrigerator in response to a torque r are given 
for an environment of two nitrogen gas reservoirs for several 
choices of the temperature T of the gases. The device works 
best when the densities of the gases, pi and p2, are equal. 
The mass of the gas particles is m = 5 x 10"^^kg, while the 
mass of the device is M = 1.66 x 10"^^kg. 



X. RESULTS AND DISCUSSION 



Eq. 70 gives us a relation between the heat flow from 
reservoir 1 to reservoir 2 and the applied torque. Earlier 
we chose the parameters of the building blocks that con- 
stitute our motor and its environment to be comparable 
with globular proteins in water, see Table|T] for details. 
Using these variables again for the refrigerating device 
we obtain a heat flow Qi^2 of 4.5 /iJ/s for every pNm of 
torque applied. This result is for T = 300 K and a shape 
that is close to optimal. Note that in the ratio (51^2/'?", 
the geometrical factor. 



M §dlrl 
TJdl^^ 



(73) 



found earlier reappears so that the comprehensive discus- 
sion therein (sections VI G and VI H ) remains applicable 



for the various shapes the parts of the refrigerator can 
adopt. 

Fig. [18] shows the dependency of the cooling power 
Qi^2/^ on the densities pi and p2 of the gas reser- 
voirs for different gas temperatures T. Here we have 
set the cooling device in the membrane separating two 
gaseous environments (the values for nitrogen gas N2 
are used: m = 5 x 10~^^kg, while for the refrigerator 
M = 1.66 X 10~^^kg and typical radius 3 nm). The 
size of the effect is determined by the ratio of the gas 
densities pi/p2- For all temperatures maximal cooling 
power is found when the densities in the two reservoirs 
are equal, pi = p2- Higher heat fluxes arise when the gas 
temperature T is higher. 

The maximal torque that can be applied and the max- 
imum obtainable temperature difference will be deter- 
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mined by increasing dissipative heat flows, which we dis- 
cuss now. 



A. Joule dissipation 

We showed the occurrence of a heat flux Qi^2 that 
takes away heat from reservoir 1 to reservoir 2. Friction 
however will cause Joule heating in both reservoirs i = 
1 , 2 by an amount 



(74) 



If we want reservoir 1 to cool down, the heat Qi^2 trans- 
ferred from reservoir 1 to 2 needs to be larger than the 
heat Qji dissipated by friction in reservoir 1, 



Qi^2 ^kBTmp2 



dirt > 1. 



Qji 1 r 

This condition poses a limit on the applied torque: 

kBTm 



r < Tiii 



P2 (p dlrf^. 



(75) 



(76) 



A greater torque would cause the heat dissipation to an- 
nihilate the cooling effect. We are able to suggest a scale- 
invariant numerical limit for the torque. Note that both 
the inertial moment / and the shape factor § dl scale 
with {R being the typical linear dimension of the re- 
frigerator) in the case that the refrigerator parts have 
homogeneous density pm- If we assume optimum opera- 
tion using the appropriate chiral shapes of the parts, as 
will be developed in detail later (section XB), we flnd 



= 1.30684. 



(77) 



This result depends only on the geometry and not on 
the dimensions of the parts. It will be lower for less 
favorable shapes and zero for a symmetric construction. 
The maximal torque than can be expressed as 



2kBT^^ X 1.30684. 



(78) 



This maximal torque is proportional to the ratio of the 
mass density p2m of gas that is heated and of the re- 
frigerator pm- It does not depend on the size of cooling 
device. For T = 300 K and a system according to Table [l| 
assuming optimum shape. 



riim = 3.92 X 10-^^Nm. 



B. Maximal net cooling 



(79) 



The cooling power of the refrigerator is proportional to 
the applied torque r ( Eg. TO] ), while the dissipative heat 
flux grows with (Eq.l74). For large r the cooling effect 



will be annihilated by dissipation and in the previous 
section we calculated a cut-off rum, at which both effects 
cancel each other. We can also calculate the torque Tmax 
that maximizes the net cooling. 



^net — Q 



1^2 



(80) 



A simple calculation leads to a maximum of Qnet 

A t - Bt'^ at Tmax Tiim/2 

(V2^^BmT)/(2/)(pip2)/((pi - 

and B = 71/7. For the optimum torque then we flnd: 



A/(2B), with A 
P2?){§dlrl)/{§dlTl) 



kBTm 



-P2 



dl rt 



(81) 



This result, like rum, is independent of the size of the 
refrigerator (assuming a homogeneous interior). It only 
depends on the density of the the environment and the 
refrigerator, the environment temperature and the spe- 
ciflc shape of the refrigerator. For an optimum shape (see 
later) and variables according to Table[l|we flnd: 



1.96 X 10"^^Nm. 



(82) 



Substituting the explicit expression for T^ax 

into Eq.[80] 

yields the maximal net heat flow out of reservoir 1: 



'net 



(fcBTm)3/2 



PipI 



{§dh 



P 



{P1+P2Y §dlrl 



(83) 



We see that when we also take into account the loss 
through friction, the refrigerator is most effective when 
the densities are the same in both reservoirs, pi = p2- 
The net cooling is higher when the device works in higher 
gas densities. The reason may seem counterintuitive: the 
heat dissipation through friction is smaller in an environ- 



ment with higher friction. Eq. 74 shows that the dissi- 
pated heat is proportional to the square of the average 
angular velocity {uo) that is obtained by applying a torque 
r, and this velocity is lower when the friction 7 is higher. 

Before we can give numerical results, we need to inves- 
tigate the role of the geometry of the refrigerator. We 
separate a scale-invariant shape factor from Eq.183) 



S = 



1 {§dlrl) 
717^ §dlrl 



(84) 



where A is the area of one part of the construction. All 
the factors that depend on the shape, the contours and 
the inertial moment, are included in S. We will begin 
straight away with optimizing this shape factor, and for 
most of the realizations of the refrigerator presented in 
this paper we assume the shape is (close to) optimum. 
The calculations however can also be done for other (less 
favorable) shapes; the elements of the calculation are 
comparable to those presented earlier, where we also dis- 
cussed three simple model shapes that can be analyzed 
analytically. A similar numerical procedure as w as use d 
to flnd the optimal shape of the motor (see section VI H ) , 



17 




FIG. 19: The effectiveness of the refrigerator is influenced 
by its (chiral) shape. S gives a size-independent measure of 
this geometrical factor. With a numerical procedure we can 
optimize this factor, by approximating the shape of the motor 
parts as piece wise linear, with n edges. For n sufficiently 
high we approach the highest obtainable factor <S, and the 
corresponding shape is represented in Fig.[20l 



mass of the refrigerator as M = 2p^A, with the 
constant density of the refrigerator parts. To find the 
maximum net heat flux the numerical factor S is multi- 
plied by 



{kBTmfl'^ 



PipI 



16 p2^A3/2 



(85) 



so that QJJ^t^ = Qo«5. Remember that our theory is two- 
dimensional. The result however depends on the abso- 
lute values of the gas densities pi and p2, contrary to 



Eq. [70] where only their relative magnitudes play a part . 
We cannot simply insert values for three-dimensional gas 
densities. Therefore we make a small technical detour 
to describe the three-dimensional expressions and justify 
for which case the two-dimensional shape optimization 
remains valid. 



C. Three-dimensional model 




In a three-dimensional analysis the essential difference 
is the description of the geometrical factors. Instead of 
a contour integral § dl there is an integral J dS over the 
surface of the refrigerator parts, while the vector r,, gets 
a new definition, see Section|V| for details. For the ratio 
of the cooling power over the applied torque we find 



Qi 



V^TTkBrnT pip2 S dS 
27 (pi+P2)^7^ 



(86) 



while the maximum net cooling power of reservoir 1 now 
becomes 



FIG. 20: The shape of one part of the refrigerating device 
that maximizes the net cooling power, which also takes into 
account the frictional dissipation. The complete device con- 
sists of this shape in one reservoir and the mirror image of 
this shape in the other. The two parts are rigidly connected 
with each other by a rotation axis, marked by a star. The 
mass distribution in the interior of the parts is assumed ho- 
mogeneous. The center of mass is represented by a dot. 



yields a value for 5, see Fig. 19 The procedure approxi- 



mates the shape of the refrigerator parts as piecewise lin- 
ear with n vertices. For sufficiently large n, S converges 
to a value slightly higher than 2. The corresponding op- 
timum shape for the cooling function is then also found, 
see Fig.[20l The shape is that of a part of the refrigerator 



in one reservoir; placing the mirrored shape in the other 
reservoir gives the optimum configuration of the refrig- 
erator. The axis of rotation is given by the z-axis. The 
numerical procedure expects the mass M of the motor to 
be homogeneously distributed, which reflects the reality 
of biological entities such as proteins that could function 
as parts of the device. 

For a homogeneous mass distribution we write the 



^ net 



8 J2 (pi+/92)2 JdSi 



(87) 



We again use a product of an external Qo and a scale- 
invariant geometrical factor tS, 



QoS. 



The geometrical factor now needs to be scaled by a factor 
proportional to R^, with R the typical linear dimension 
of the shape. For this we use V^/^, with V the volume 
of a refrigerator part: 



1 M^ iJdSr 
^27372 JdSi 



(89) 



The shape factor can then be solved analytically or nu- 
merically, producing size-independent results. The re- 
maining factor 



/tt {kBTm) 



3/2 



P1P2 



M2 



{pi + P2)' 



_y2/3 



27r (kBTmf/'^ 



PipI 



16 p^FV3 {pi^p^Y 



(90) 
(91) 
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For a globular pro- 
3nm, we find a value of 



2.992 X 10~^^kg for an aquatic environment, and pm = 
1380 kg m~^, typical for proteins, Q^q^ = QoS can 
be expressed as a function of only the radius R of one 
refrigerator part, as in Fig. 21 
tein of typical dimension R 
Q^&t = 3.5 X 10"^^J/s. Note the strong size depen- 
dence: a refrigerator of one nm, would yield a cooling 
power of about 2.8 x IQ-^^J/s. For comparison, it takes 
about 2.2 X 10~^^J to cool down a 1 /am small cell one 
Kelvin, which could be accomplished by one refrigerator 
of 1 nm radius in one minute. 



FIG. 21: The net cooling power of the Brownian refrigera- 
tor (the dissipative loss by friction is deducted from the heat 
flux that cools down one compartment) depends strongly on 
the size of the cooling device. Values for a device operat- 
ing between two aquatic reservoirs (m = 2.992 x 10~^^kg, 
= P2 = 3.343 X lO^^m"^) at temperature T = 300K are 
presented here. The size-dependency of the mass of the rotat- 
ing device M is included, while its density is that of a typical 
protein, pm = 1380 kg m~^. 



then shows again the role of the various parameters. In 
the second expression for Qo we substituted M = 2pmV 
for the homogeneous case. 

We now argue that we can recuperate the two- 
dimensional optimization of the shape factor S and corre- 
sponding two-dimensional shape. We propose a prismatic 
structure for each of the refrigerator parts, defined by two 
equal, flat (two-dimensional) surfaces separated by a dis- 
tance as shown in the initial model of the device in 
The rotation axis is perpendicular to the two 



Fig. 17 



surfaces. For this configuration it is easy to show that 
the surface integrals can be written as the product of the 
distance H and the contour integral of the top or bottom 
surface, exactly as in the two-dimensional description: 



H (bdl K 



dSrt =H 



d^ ^%2D- 



(92) 
(93) 



(94) 



Also, for a homogeneous mass distribution, 

M _ V _ HA 
T ~ jr\dV ~ Hjr^dA' 

where we take the integral over the volume V of the 
distance to the rotation axis, and find it is equivalent to 
taking the surface integral over the top or bottom sur- 
face A of the prism times the thickness H. In conclu- 
sion we recover the two-dimensional shapefactor (Eq.[84|) 
by inserting the corresponding volume V = (i^V^)^ 
The approach is then to use the three- 



Eq. 89 



dimensional expression for Qq (Eq. 91), and the numeri- 
cal results for the two-dimensional case, as in Fig. [19] 

Doing this for a temperature (T = 300K), reser- 
voir densities {pi = p2 = 3.343 x lO^^m"^) and m = 



D. Thermal conductivity 

As mentioned earlier, the L22 Onsager coefficient can 
be related to heat conductivity between the reservoirs. 
In the linear response model, the heat conducted from 
one reservoir to the other can be quantified: 



Qcond = ^^/C^AT. 



(95) 



Here 71 and 72 are the friction coefficients of the sepa- 
rate parts of the refrigerator, and 7 = 71 + 72 represents 
the total friction coefficient. The conductive heat flow is 
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FIG. 22: The temperature gradient AT for which the conduc- 
tive heat flow cancels the cooling heat flow. For the device 
to cool down reservoir 1 (at the expense of reservoir 2), the 
temperature difference must remain under this limit. The 
different curves are for examples of the mass ratios between 
colliding particle (m) and device M. A small M (or large m) 
is beneficial. The sustainable temperature gradient is highest 
for high gas density p2 in reservoir 2 (relative to pi). For 
high gas density pi in reservoir 1 the obtainable temperature 
gradient becomes very small. The graph is for absolute tem- 
perature differences around T = 300 K. 

proportional to the temperature gradient AT and goes 
from the warm to the cold reservoir. Therefore it is di- 
rected against the cooling power of the device and for a 
temperature difference larger than ATum the net cooling 
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effect will vanish. The condition Qcond 

< QJfet'' l^a^is to For a bar-shaped motor, a = or a = 90°, Eq. dA3| ef- 



AT irm 



P2 



{§dlrl 



v2 • 



(96) 



Note that the sustainable temperature gradient is pro- 
portional to m/M. In Fig. 22 we show ATum for varying 



m/M as a function of the density ratio pi/p2- For our 
choice of the direction of the cooling effect (reservoir 1 
cools down, reservoir 2 heats up), a higher relative den- 
sity in reservoir 2 allows a larger ATum. In the limit 
of pi being much smaller than p2, ATum/T approaches 
{7r/S){m/I) [§ dlrf^)^ / [§ dlr^)^ Earlier we found equal 
densities, pi = P2, to correspond with maximum net cool- 
ing power, for which the maximum sustainable tempera- 
ture gradient ATum is half this value. 

For the previously used example of globular proteins 
(see Table[l|for the parameters), ATum = 4.3 mK at T = 
300 K and equal reservoir densities. 

In all examples given, the shape is presumed to be 
(close to) optimal. 



APPENDIX A: FULL EXPRESSIONS FOR 
MODEL MOTORS 

For Motor 1 a straightforward calculation leads to 
R 



dl r,^ = — f cos 4a 
" 72 V 

+ 2i? ((i? + ^y) cos 3a - (i? + 6x) sin 3a) 

+ 12 ((?/ — x) (i? + 3(x + y)) cos 2a — Qxy sin 2a) 

+ 6((i?^ + 2^i?+ 12^2)cosa 

+ (i?^ + 2xi?+ 12x^)sina) 

+ IR^ + I2{x + y)R + 36(x^ + ^^) ) , (Al) 



while I/M = R'^/IS + + y'^. The friction coefficient is 
then directly obtained through Eq. (29). One can note 



that this coefficient is always nonzero, unless R = 0. Also 
for a = 0, when the motor has the shape of two connected 
bars, Eq. ( |A1[ ) is nonzero. Choosing a = 0, x = 0, and 
y = i^/3, we retrieve the expression for a rod of length 
R rotating about one end. 



dlr'f^ = 2R^/3. 



(A2) 



The sense of rotation is geometrically determined by 
Rsin 2a 



dl 



432 



R"^ sin 4a 



ISR^ (x cos 3a -\- y sin 3a) 
10SR{2xy cos 2a - (x^ - y'^) sin 2a) 
M{x{R^ -4x^ + 12?/^) cos a 

y{R^ + 12x^ - 4^^) sin a) ) . (A3) 



fectively returns zero, resulting is a zero average rotation 
for the motor. The shape needs to have chiral asymme- 
try for the motor to function. This observation becomes 
more apparent when we locate the rotation axis in the 
center of mass, (x^y) = (0,0). Then Eq. (A3) simply 
becomes 

= - ^ sin(2a) sin(4a) , ( A4) 

which equals zero also when a = 45°. For this angle the 
shape contains an extra symmetry axis, lifting the chiral 
symmetry. To find a nonzero rotation, we need to change 
the shape to a chiral asymmetric one. We see an opposite 
rotation sense when moving across the symmetric posi- 
tion. In particular for this choice of the position of the 
rotation axis, we note that {u)[7r/2 — a] = — (cj)[a], where 
for this triangular shape, a and 7r/2 — a correspond to 
chiral enantiomers. This is because §dlrf^^ ^/I^ and 
vq are even under aa^7r/2 — a transformation, while 
§ dl rf^ is uneven. For other locations of the rotation axis, 
we see a similar behavior. In Fig. |9] we show the depen- 
dency of the angular velocity (cj) on the exact place of the 
axis (which in our calculations can be outside the motor 
body). For chiral shapes, the axis of rotation that yields 
zero rotation is no longer on a symmetry axis. These 
locations are generally given by three curves in the two- 
dimensional plane. Nonzero rotation is found in the six 
zones divided by the zero curves and the sense of rotation 
is opposite in neighboring zones. One point of maximal 
rotation speed is found per zone, actual and maximal 
values of {uo) will be discussed later. To conclude our 
discussion of Motor 1, a few expressions for fixed values 
are given. For a = 7r/4, 



dl r,^ 



R^ 
36 



(3 + 2^2^ R^ + 6{x + y)R 
18 ((1 + 72) x^ - 2yx + (1 + 72) y^^y (A5) 



R 

dlrl = Y^{x-y) 

X {q{x + y)R + 6^2 (x^ + Ayx + y'^) - V2R^^ , (A6) 
and for a = 7r/6 

dlr^ ^ 



= — ((5 + 2V3) i?2 + 4 (3 + v^) yR 
36x^ + 12y (^(3 + 2Vtj y - 2\/3a;) j , (A7) 



- mV3{x -y){x^ y)R - 27y {R^ + 12x^ - 4?/^) 

+ 27V3x {R^ - 4x^ + 1 V)^ . (A8) 



20 



Similarly for Motor 2 we find, again choosing the origin 
of the coordinate system (x^y) in the center of mass of 
the motor, 



2,R [ 

dl r,^ = — I + 3Ry cos 3a + 9y^ + 9x^ sin a 
9 V 

sin^ a [R^ (3 sin a — 2 cos 2a) + 9x^-9^^) ), (A9) 



while I/M = R^{2 - cos2a)/18 + + y^. Again 
^ dl rf^ = sin 2a (^R{3y + R cos a) (1 — 2 cos 2a) 

+ 3(x^ -3?/^)cosa^ (AlO) 

determines the sense of rotation and reflects the chiral 
symmetry of the motor. Specifically, putting the rotation 
axis on the symmetry axis of the motor, x = 0, results in 
a zero average rotation. Putting the axis on symmetrical 
location across the symmetry axis yields the same rota- 
tion speed, but in opposite sense, (uj) [— x, y] = — {uj) [x, y] . 
Again, reducing the motor to a bar shaped object, a = 
or a = 7r/2, a preferred sense of rotation can no longer 
be determined and {uo) — 0. Fig. 10 shows (cc;) for general 
positions of the rotation axis, for different configurations 
of the motor shape. For certain configurations the expres- 
sions simplify to reflect a higher symmetry. For a = 7r/6 
for example the shape is that of an equilateral triangle 
with side 2R^ with 



dl rt 



1 



(All) 



and 



j dl rl = ^Rx {x^ - 3y^) , (A12) 

so that the full expression for the rotation speed becomes 
36x (x^ - 3y'^) 



R^ + 18 (x2 + y^) R^ + 72 (x2 + y^y 



(A13) 



The rotation speed will be zero, when the rotation axis 
is on any of the three symmetry axes of the equilateral 
triangle, x = ±V3y besides x = 0, while the extrema 
are located on the vertices of a regular hexagon of side 



(3 + v33)/6i?. Finally, for Motor 2, the expressions 
for a = 7r/4 are 

jdlrl = ^(^(4.^ 3V2) R^ - 6V2yR 

+ 18((l + y2)x' + ^') j, (A14) 



dl rt 



Rx 



(y2R^ + 6yR + 3V2 (a 



(A15) 



For Motor 3 the origin of the coordinate system (x, y) 
is chosen in the center of the disk sector, and the expres- 
sions for a general shape (a), size (i?), and location (x, y) 
with respect to the rotation axis are 

^ dl rl = ^ (^2R^ + sin^ a + 3 (x^ + y'^) a 
+ 6y cos a{y cos a — R) -\- 3{x — y){x -\- y) cos a sin a^ , 

(AI6) 

even for transformation, and 

dl rt = ^ (x^ - 3y'^) sin 3a 

— 3R{R — 3y cosa) sinaj ^ (All) 



odd for a X ^ X transformation, so that the sense of 
rotation is inverted by putting the axis of rotation at its 
mirror location with respect to the symmetry axis, x = 0, 
of Motor 3. This is because the inertial moment, 



R'^ ARysina ,2,2 



(A18) 



is also even for the x — x transformation. For certain 
cases the expressions become simpler, such as for a semi- 
disk, a = 7r/2: 

^ dl rl = I (4i?2 + 3(4 + 7r)x^ + 3 V) , (A19) 



and 



dli 



2Rx 



(3i?^ +x^ -3^^) , (A20) 



and for a = 7r/3: 



dl rt 



R 
12 



(^R^ - 12yR + (18 + 3^3 + Att^ x^ 

+ (6-3y3 + 47r)^2j, (A21) 



and 



dirt = ^R^x{3y-2R). 



(A22) 



Fig. 11 shows (cj) for these and other configurations as 
a function of the position of the rotation axis. Again 
extrema are obtained in the x, ?/-plane at finite values, 
while three curves of zero rotation mark areas of opposite 
rotation. 



[1] R Reimann, Phys. Rep. 361, 57 (2002). 
[2] R.D. Astumian, Science 276, 917 (1997). 
[3] F. Jiilicher, A. Ajdari, and J. Prost, Rev. Mod. Phys. 69, 
1269 (1997). 

[4] S. Leibler, Nature (London) 370, 412 (1994). 

[5] M. von Smoluchowski, Phys. Z. 13, 1069 (1912). 

[6] R. Landauer, J. Stat. Phys. 53, 233 (1988). 

[7] C. Van den Broeck, R. Kawai, and P. Meurs, Phys. Rev. 
Lett. 93, 090601 (2004). 

[8] P. Meurs, C. Van den Broeck, and A. Garcia, 
Phys. Rev. E 70, 051109 (2004). 

[9] C. Van den Broeck, P. Meurs, and R. Kawai, New J. 
Phys. 7, 10 (2005). 
[10] P. Meurs and C. Van den Broeck, J. Phys.: Con- 
dens. Matter 17, S3673 (2005). 
[11] J. M. R. Parrondo and P. Espagnol, Am. J. Phys. 64, 
1125 (1996). 

[12] K. Sekimoto, Progr. Theor. Phys. SuppL 130, 17 (1998). 
[13] C. Jarzynski and O. Mazonka, Phys. Rev. E 59, 6448 
(1999). 

[14] C. Van den Broeck, Adv. Chem. Phys. 135, 189 (2007). 
[15] R. D. Astumian, Phys. Rev. E 76, 020102 (2007). 
[16] C. Van den Broeck and R. Kawai, Phys. Rev. Lett. 96, 
210601 (2006). 

[17] N. Nakagawa and T. S. Komatsu, Europhys. Lett. 75, 22 
(2006). 

[18] J. P. Pekola and F. W. J. Hekking, Phys. Rev. Lett. 98, 
210604 (2007). 

[19] L. Onsager, Phys. Rev. 37, 405 (1931); ibid. 38, 2265 
(1931). 

[20] H. B. Cahen, Thermodynamics and an Introduction to 
Thermo statistics (Wiley, New York, 1985). 

[21] H. Noji, R. Yasuda, M. Yoshida, and K. Kinosita, Na- 
ture 386, 299 (1997). 

[22] R. Yasuda, H. Noji, K. Kinosita, Jr, and M. Yoshida, 
Ceh 93, 1117 (1998). 

[23] H. Itoh et al.. Nature 427, 465 (2004). 

[24] R.P. Feynman, R.B. Leighton, and M. Sands, The Feyn- 
man Lectures on Physics I (Addison- Wesley, Reading, 
MA, 1963). 

[25] N.G. van Kampen, Stochastic Processes in Physics and 
Chemistry (North-Holland, Amsterdam, 1981). 

[26] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, 
Series, and Products (Academic, New York, 1980). 

[27] M. van den Broek and C. Van den Broeck, 
Phys. Rev. Lett., to appear. 



